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We present a theory describing spiral magnets with Dzyaloshinskii-Moriya interaction (DMI) 
subject to bond disorder at small concentration c of defects. It is assumed that both DMI and 
exchange coupling are changed on imperfect bonds. Qualitatively the same physical picture is 
obtained in two models which are considered in detail: B20 cubic helimagnets and layered magnets 
in which DMI leads to a long-period spiral ordering perpendicular to layers. We find that the 
distortion of the spiral magnetic ordering around a single imperfect bond is long-range: values of 
additional turns of spins decay with the distance r to the defect as 1/r^ being governed by the 
Poisson’s equation for electric dipole. At finite concentration of randomly distributed imperfect 
bonds, we calculate correction to the spiral vector. We show that this correction can change the 
sign of spin chirality even at c <C 1 if defects are strong enough. It is demonstrated that impurities 
lead to a diffuse elastic neutron scattering which has power-law singularities at magnetic Bragg peaks 
positions. Then, each Bragg peak acquires power-law decaying tails. Corrections are calculated to 
the magnon energy and to its damping caused by scattering on impurities. 

PACS numbers: 75.10.Jm, 75.10.Nr, 75.30.-m, 75.30.Ds 


I. INTRODUCTION 

In crystals without center of inversion, Dzyaloshinskii-Moriya interaction (DMI) is caused by an antisymmetric 
spin-orbit interactioiiiii^ The competition of the symmetric ferromagnetic (FM) or antiferromagnetic (AF) exchange 
interaction and DMI can result in a spiral magnetic structured Although a long time has passed since the spiral 
ordering was observed for the first time, helimagnets with DMI still attract a lot of attention. This interest is 
stimulated by discovery of rich phase diagrams and exotic spin structures caused by DMI which arise under certain 
conditions. Phases with such topological states as chiral soliton lattices in layered helimagnets (e.g., in Cri/ 3 NbS 2 )^ 
and skyrmion lattices in B20 cubic chiral magnets (e.g., in MnSi)^ are widely discussed now. These materials are 
attractive not only from a fundamental but also from a technological point of view owing to their potential applications 
in spintronic devices. 

Mixed B20 spiral compounds have been considered experimentally recently^ It is shown in Refd that the modulus 
of the spiral vector q in Mni_ 2 ;Fea;Ge depends on dopant concentration x and the magnetic chirality changes its sign 
(and q goes through zero) at a: « 0.75. This observation is quite expected because MnGe and FeGe are B20 cubic 
helimagnets with opposite signs of the spin chirality. Evidently, such a behavior is a consequence of the fact that the 
exchange interaction and DMI change around dopant ions which can be considered as defects at a: <C I or a: ~ 1. 
These experimental results are interpreted phenomenologically by renormalization of constants in the Hamiltonian 
describing the pure translationally invariant B20 magnets. Then, a more detailed theoretical description of mixed 
spiral materials is required. 

Motivated by this experimental activity, we address in the present paper the problem of spiral magnets with DMI 
subject to bond disorder at small concentration c of defects. We assume that both exchange interaction and DMI are 
changed on imperfect bonds. Two models are considered in detail: i) spiral magnets on a simple cubic lattice with 
FM exchange coupling and small DMI between nearest-neighbor spins, where the DMI vector is directed along the 
line connected couple of spins, and ii) layered magnets with small DMI which acts between nearest-neighbor spins 
from different layers and which vector is directed along z (chiral) axis perpendicular to layers (see Fig. [T|). The most 
famous and the most studied compounds described by the model of the first type is probably MnSi and those of the 
second type are Gri/ 3 NbS 2 and CsGuGl 3 . At zero magnetic field and small temperature, DMI leads to long-period 
helix structures in these materials along one of the space cubic diagonals and along z axis, correspondinglyii"— 

We obtain qualitatively the same physical picture in both models. The one-impurity problem is addressed first. We 
show that the perturbation of the spiral ordering around the defect bond (i.e., values of additional turns of spins due 
to the defect) is described by the Poisson’s equation for electric dipolei^ Then, the magnetic ordering disturbance 
made by one impurity is long-range: values of additional turns of spins decay with the distance r to the defect as 1/r^. 
This finding can be easily extended to the corresponding models on lattices with space dimensions d > 2, the result 
being The spin texture around a ferromagnetic bond observed in two-dimensional collinear AFs follows the 

same law (see Refs j^°d^ and references therein). It has been found recently that the distortion of magnetic ordering 
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FIG. 1: (Color online.) Two types of spiral magnets with DMI considered in the present paper, (a) Layered spiral magnet ([T]) 
with tetragonal lattice in which DMI acts only between nearest-neighbor spins from neighboring xy planes (DMI vector D is 
depicted). Exchange coupling constants between neighboring spins inside xy planes (Ji) and along « axis (Jo) are also shown. 
The long-period helix propagates along z axis. Similar models with hexagonal xy planes describe Cri/ 3 NbS 2 and CsCuCls. 
The imperfect bond is shown by dashed line, (b) Cristal structure of MnSi that is probably the most famous representative of 
spiral cubic B20 magnets. The helix can propagate along any space diagonal of the cube. 


around defects decays exponentially in collinear AFs in magnetic fieldi^rJd while the spin texture around the vacancy 
in triangular AF decay as 

At finite defects concentration c <C 1, spiral magnets we discuss are equivalent to a dielectric with randomly 
distributed electric dipoles which lead to a finite average “polarization” of a unit volume proportional to c. This 
“electrical polarization” corresponds to a correction 5q do the modulus of the spiral vector q. 

Our analysis of the elastic neutron scattering cross section predicts magnetic Bragg peaks (satellites) on momenta 
transfer Q = ±(q -h i5q) + r, where t is a reciprocal lattice vector. Besides, we obtain a diffuse scattering. Quite 
unexpectedly for diffuse scattering caused by disorder, its cross section has power-law singularities at positions of 
magnetic Bragg peaks. This feature is attributed to the long-range character of the perturbation made by defect 
bonds. Thus, impurities result in the shift by (5q of the magnetic Bragg peaks positions and in power-law decaying 
tails of each peak. 

We calculate also magnon spectrum renormalization due to the scattering on defects in the first order in c. These 
calculations are performed in the layered helimagnets only for FM exchange coupling constants. 

The rest of the present paper is organized as follows. Secs. [TTl and mi which have similar structures are devoted 
to layered helimagnets and to B20 cubic spiral magnets with DMI, respectively. In subsections [TIlA and IIIIA . we 
consider Hamiltonians of pure systems using the conventional Holstein-Primakoff spin transformation. We discuss in 
subsections and the perturbation of the magnetic ordering around one imperfect bond and consider small 
concentration of such bonds. Then, we present our results for elastic neutron scattering cross-section in systems 
with bond disorder (Secs. Hit ) and HTTtD) . In Secs. HID and HTIt). magnon spectrum renormalization is considered. 
All calculations in Secs. HK-Et) are carried out for layered helimagnets with FM exchange interaction. We show 
in Sec. HI El that these results (except for the spectrum renormalization) are applicable after simple modifications to 
many other layered helimagnets with bond disorder. Sec. IIVI contains the summary and the conclusion. One appendix 
is added with some details of the magnon spectrum calculation. 


II. LAYERED SPIRAL MAGNETS WITH DMI 
A. Pure system 

In this section we consider a magnet containing FM xy planes with a simple square lattice and the exchange coupling 
between neighboring spins only. Planes are stacked along z axis. We take into account the exchange coupling and 
DMI between neighboring spins from neighboring planes. The DMI vector D = De^ is the same for all bonds along 
z axis, where = (0, 0,1) is the unit vector directed along z axis and we assume for simplicity that the distance 
between all neighboring sites is equal to unity (see Fig. [Ija)). The Hamiltonian of this system has the form 

— To ^ ^ J\ ^ ^ ^ ^ D • X , (1) 

in (ij)n in 
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where Jo, Ji > 0, Jq, Ji 3> -D, Si„ is an operator of the spin sitting at i-th site of n-th plane, and {ij)n denote nearest 
neighbor sites in the n-th plane. The last term in Eq. o containing antisymmetric combinations of spins x Si„_|_i] 
can be eliminated by applying the rotation about z axis by a pitch 

Sfn = Sf^cosnq-Sf^ sin nq, 

Sfn = S^^smnq +cos nq, (2) 

QZ qz' 

^in ^in’ 

The value of q is chosen so that the antisymmetric spin combinations disappear in the Hamiltonian. Simple calculation 
give 


tang = — <C 1. 
Jo 

After transformation @, Hamiltonian ([1]) obeys the following form: 


(3) 


n = - 


[jcsi 




+ + 


— Jl 


E 


S' S' 

jn") 


(4) 


in 


{ij)n 


where Jq = Jo\/l + {D/ JqY and S'„ = S'f„, S'f„). Thus, the initial Hamiltonian ([T]) of the system with the 

spiral spin ordering described by the vector q = (0,0, g), where g is given by Eq. (jS]), is equivalent to a EM described 
by Hamiltonian As Jo > Jo, xy plane is the easy one in EM (|3]). Then, DMI forces spins to lie in the plane 
perpendicular to D. 

Eor further consideration of Hamiltonian (|4]) , we use the Holstein-Primakoff spin representation 


QX 


qy 

^in 





atnOin 

AS 


^tn ^in \ 
AS ) ’ 



atnOju 

AS 


+2 \ 
45' ) 


(5) 


After simple calculations, one obtains that there are no terms in the Hamiltonian which are linear in Bose operators 
and which contain products of three Bose operators. Terms containing products of two operators of creation and 
annihilation have the form 


1~L2 — SJq ‘2SJ\ ^ ^ ^ 


n^^Clin 




( 6 ) 


where we omit terms of the second order in D/ Jq <C 1. 


B. Perturbation of the magnetic ordering by defects 


Let us discuss a defect bond with DM vector D' = (0, 0, D') ^ D and Jq ^ Jq between spins at sites 00 and 01 (see 
Fig. UKa)). The following additional terms arise in Hamiltonian O: 


V = Vcim H“ Vea; = Udm^z * [SqO ^ Sgi] '^^eaiSoo ‘ Sgi, 
Udm — D Z?, 

'^ex ~ Jq- 

One obtains for the perturbation of Hamiltonian Q from Eq. © using Eqs. m and dSI) 


Vrfm — SudjYi^j ^ S Ctoi ^00 ^Oo), 


Tea: — Site 


2 

SD_ 

2Tq 


-I- aoi - a(}o - ooo) -b Suex{a^iaoi + aQQaoo - aj^aoo - a^JpOoi), 


(7) 

( 8 ) 

(9) 


( 10 ) 

( 11 ) 
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where we take into account only linear and bilinear terms in Bose-operators which are of the zeroth and of the first 
orders in DMI (in particular, we put cosq = 1 and sing = D/Jq). 

Terms in Eqs. (HOD and dm linear in Bose-operators signify a distortion of the FM ordering around the imperfect 
bond. To eliminate the linear terms in the Hamiltonian, one has to make the shift 

^in — ^in T i (1^) 

“in “ ^in ' Pmf> j 

where and ipin are constants (the “condensate density” and the “phase”, respectively) which describe perturbation 
of the spin ordering due to the defect. As the easy-plane anisotropy in the Hamiltonian forces all spins to lie within 
the xy plane and we do not consider AF coupling on the defect bond (i.e., Jg > 0), we put (pm = 0 in the following 
to eliminate the magnetization component perpendicular to the easy axis. As it is seen from Fq. ([S]) and illustrated 
by inset in Fig. [UJa), a real pm 7 ^ 0 describes a rotation of the magnetization at site in within xy plane. To restrict 
ourselves to terms of leading powers in pin in subsequent calculations, we assume that \pin\ '/S. Then, the rotation 
angle is equal approximately to pin\/^2jS (because ~ V^pin and = S — pf^K. S). 

Bilinear part of the Hamiltonian ([5]) acquires the following form after shift (fT^ : 

H 2 = (13) 

= SJg (2pin - Pin+I - Pin-i) + 25'Ji 6+ (p^ - Pjn) + h.c., (14) 

in {ij)n 

where h.c. denote the Hermitian conjugated terms, does not contain Bose operators, and is obtained from 
Eq. ([ 6 |) by the replacement of operators a by operators b. One has to dispose of linear in bin and bf^ terms in the 
Hamiltonian by choosing proper As usual, a minimum of the classical energy (i.e., the part of the Hamiltonian 
not containing Bose-operators) is realized at those pm which cancel the linear terms in the Hamiltonian. Let us find 
such Pin in two steps: we consider first Vdm only assuming that Uex = 0 and then we take into account both Vdm and 
Vex given by Eqs. (US and (EH), respectively. 


1. Defects in DMI only (uex = 0 ) 


We start with the one-impurity problem and then we consider a finite concentration of defects. One has from 
Eq. (fTUl) after shift (IT^ 

Fs 

Vdm = UdmSy —(2poO ~ 2poi + ■)" ^00 ~ ^01 “ ^Ol)- (15) 

In order linear terms die out in the Hamiltonian, the following equations should hold for every site in which follows 
from Eqs. El and El 

dj ^ ' {pin Pjn) Jo{‘^Pin Pin — 1 Pin+l) — tldmy ~^{^in,00 ^m,0l), VzTT, (15) 


where j enumerates nearest neighbors of i-th site in n-th plane and <5 is the Kronecker delta. It is well known that 
the second derivative of a function f{x) can be written as 


f{x + h) + f{x - h)- 2f{x) 


dx"^ 




(17) 


with a good precision if f{x) does not change considerably at a distance of h. Thus, Eqs. (ITdl) can be represented in 
the differential form in the continuum limit as follows: 

Da 

(18) 


fd‘^p{r) d^p{r)\ d‘^p{r) 


where (5(r) and d(r—rg) are delta-functions defining positions of two spins involved in the defect bond and rg = (0,0,1). 
One expects that the solution of Eq. (IT^ describes well the solution of Eqs. (fTHll not very close to the imperfect bond, 
in which region p(r) changes rapidly. After rescaling in xy planes 




Jo 


X = \ —X and y = \ -rU^ 


Ji 


Ji' 


(19) 
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FIG. 2: (Color online.) (b) xz plane containing the defect bond which is shown by bold line. Lattice sites are drawn by blue 
circles. Contour plot is also shown at Ji = Jo/4 of the function p(r) given by Eq. (1211) and divided by Q/Jq. (a) and (c) 
Plots of p(r) along dashed lines depicted in panel (b). Squares represent the result of numerical solution of Eqs. (1161) for the 
cluster with 24 x 24 x 24 sites. It is seen that the analytical result (pT1 ) -(|2 ^ starts working right from sites adjacent to the 
defect bond. Inset in panel (a) illustrates the meaning of the “condensate density” pi„ (as well as its counterpart p(r) in the 
continuum limit) in our consideration. Appearance of pin yf 0 at site in signifies a rotation of the magnetic moment at that 
site by angle pin a/2 /S in the xy plane. 


Eq. (UHl) turns into the Poisson’s equation 


Ap(r) 


'^dm 

Jl 


|(5(r) -(5(r-ro)), 


( 20 ) 


where A = d'^jdx'^ + jdy^ + d^jdz"^. Eq. (|2ni) describes electrostatic field of a dipole 


p(r) 

Q 

d 


Q 

( 1 

47r Ji 

llf-rol 



'^dm \ 

2 ’ 

1 Q 


47r Ji 



( 21 ) 

( 22 ) 

(23) 


where f = (J, y, z) and d is the dipole moment. Then, the magnetic ordering distortion produced by one defect bond 
is long-range: it decays with the distance r as 1/r^. We observe by numerical solution of Eqs. (HU) that the result 
dUD starts working well right from sites neighboring to the defect bond in a broad range of parameters (see Fig. m- 
The correspondence between the electrostatic picture and the spiral magnet ([T]) is illustrated by Fig. [H 
Distortion of the EM ordering in the spin system ([4]) with a finite concentration c <C 1 of such randomly distributed 
defects is described by the electric field from a set of randomly distributed dipoles having the same dipole moment d 
given by Eq. (ESI). Averaging over the system volume, one obtains for the “electric polarization” 


P = (24) 

Vo 

where vq = Jq/Ji is the unit cell volume after rescaling (|T^ . The field p(r) inside the uniformly polarized system is 
given by the equation 

Vp(r) = 47rP 

which has the following explicit form in our case: 

&p{r) _ c Udm _ Udm 1^ 

dz Do Ji V 2 Jo V 2 

The solution of Eq. (1261) gives an averaging solution of our problem which has the form 


(25) 

(26) 


p(r) 



(27) 
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FIG. 3: (Color online.) Illustration of the correspondence between the electrostatic picture proposed in the text and the spiral 
magnet ©• Horizontal dashed line represents the plane perpendicular to the defect bond and to the dipole moment (I23II . The 
condensate density p(r) given by Eq. (1211) (the field of the dipole) has opposite signs above and below this plane. This signifies 
that spins lying above and below this plane acquire additional turns in opposite directions which are depicted by blue arrows 
and which values are governed by |p(r)| as it is explained in the text (see also inset in Fig. [2ja)). 


where we omit a constant corresponding to a rotation of all spins in the system by the same angle. Eq. (EZl) corresponds 
to the following correction to pitch (jS]): 


Sq 


qV 


St 


QXf 

^in—1 


= C- 


'^dm 


< 1 . 


(28) 


It should be noted that the requirement \pin\ yfS is essentially important for calculations leading to Eq. (EZD 
whereas Eq. (EZD contradicts it. This discrepancy can be easily removed by applying rotation m by the pitch q + Sq 
rather than by g, where q and Sq are given by Eqs. ([3]) and (I28L respectively. Carrying out again the corresponding 
calculations, we obtain, in particular, Eq. (ESI) with P = 0 which solution is p(r) = 0. Then, we conclude again (not 
violating the requirement \pin\ "C V^) that defects lead on average to the correction (1281) to the spiral pitch ([3]). 

After rotation (El) by the pitch q + Sq^ one obtains for the condensate densities not very close to impurities 


= E- 

1=1 


JV. d • (r„ - R,) 




- dTTC^d • R„, 
Jo 


(29) 


where d is given by Eq. (l23l) and j enumerates Nd defect bonds in the system. The first term in Eq. (l29l) is the field 
from all dipoles in the lattice and the second one arises due to the additional turn by Sq. Naturally, averaging of 
Eq. (1331) over the whole system gives zero. 


S. Defects both in the exchange interaction and DMI 


Taking into account also the imperfection of the exchange interaction m on the defect bond, we obtain from 
Eqs. ( 0 , (nni), and EH) for the part of V which is linear in Bose operators 


(30) 


— S — Uex-J-^ ~ POl)^ (^00 + ^00 “ ^01 “ ^Ol) ■ 


The counterpart of Eq. (HE) has the form in this case 


jPin Pjn) S' J{i{f2.piji Pin —I Pin^i^l) — (*^m,00 ^m,0l) 


D 

'^ex “7 '^dm 

Jo 


UexiPOl — Poo) 


\/in. (31) 


3 L - ' J 

These equations are more complicated than Eqs. (1161) because one cannot solve them directly in the continuum limit. 
As it is pointed out above, the solution in the continuum limit does not describe the solution of the initial equations 
near the defect bond. On the other hand, the solution of Eqs. (I31|) in the continuum limit is essentially determined 
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by condensate densities at sites involved in the defect bond (because the right-hand side depends on poi and poo)- 
Then, we use the following self-consistent scheme to solve Eqs. (EH). First, we put 

Poi — Poo = o. (32) 

in the right-hand side of Eqs. (1311) and treat a as an unknown constant. As a result one returns to the problem 
considered in the previous subsection which solution is given by Eq. m, where now 


Q — ^2 ^ex j ^ 

Second, we consider two equations EH for n = 0 and n = 1 

4 Jl(poO — Pio) + <70(2^00 — POl — PO-l) = —Uex{p00 — P0l) — 


4Jl(poi “ Pll) + 7o(2poi — Poo ~ P 02 ) — Uexipoo — Poi) + \/~ ( "^dm — Uex-J- ) , 


'^dm '^ex j J ) 


D 


( 33 ) 

( 34 ) 

( 35 ) 


where pio and pn are condensate densities at sites neighboring to spins involved in the defect bond and lying in 
planes with n = 0 and n = 1, respectively (i.e., we use the system symmetry). One obtains by subtracting Eq. (IMl) 
from eh 


Ji( 4 a - 8pii) -I- Jo( 3 a - 2 po 2 ) = - 2 aUex + (^Udm - '^ex^^ ’ 


(36) 


where we use that one can turn the coordinate system to fulfill relations po 2 = —Po-i and pio = —pii. Using our 
finding that the result (1211) obtained in the continuum limit starts working well right from sites neighboring to the 
defect bond in a broad range of parameters, we derive po 2 and pn from Eqs. EH and (1551) . Then, Eq. (|55)) turns into 
an equation for a which solution is given by 


a = 



Udm Uex 


Jo J 


2 


3To 4Ti -)- u 




One has from Eqs. (1551) and EH 


Q = 



D 

'^dm '^ex “7” 

do 


3t/o 4(/i ~h Uex 


3Jo + 4(7i 

_|_ dp \ 2 l f / 

47 rJi 7 T \y Jq Y Ji + Jo J 


( 37 ) 


( 38 ) 


It should be noted that Eq. EH can give an infinitely large result at Ji < O.Uq and — 1 < Uex < 0. This signifies that 
more equations (1511) should be considered in addition to Eqs. (l34)) and (1351) to find poo and poi. The corresponding 
analysis is out of the scope of the present paper. Our numerical solutions of Eqs. EH on finite clusters show that 
Eqs. (15X1) and (1551) work well beyond the region Ji < O.lJolJ—1 < Uex < 0 not very close to the defect bond. 

As it is done above, we derive for the correction to the spiral pitch (cf. Eq. (1551) 1 


6q 


D 

'd'dm 

do 


3 Jo -|- 4 Ji -|- Uex 


3 + 4 Ji/Jo _ 

_|_ Jo \ i 

47 rJi Jp V + ) 


(39) 


Values of pm are determined by Eq. ([29|) . where d and Q are given by Eqs. ([23)) and (155|) . respectively. 

It is interesting to note that the influence of imperfections in the DMI and in the exchange interaction on the spiral 
ordering weaken each other substantially at Udm ~ UexD/J^ (see Eqs. (|38)) and (1591) '). 


C. How defect bonds seen in elastic neutron scattering experiments 

The cross-section of elastic neutron scattering is given by^^ 

m,jm X:V 


( 40 ) 





















where Q is the momentum transfer, Q = Q/Q, = x,y, (...) denotes an average over quantum and thermal 
fluctuations, 

i^m) ~ ^ COS nq' — V^pin sin nq' — cos nq', (41) 

{Sfri) ~ Ssuvnq' + Pin cos nq' - cos nq', (42) 

q' = q + 5q, (43) 

Pin are given by Eq. q and 5q are given by Eqs. (jS]) and (1551) . respectively, and we omit terms of orders higher 

than the second power of p. Terms in Eqs. (HD) and (H5|) not containing p lead to the well known result for a spiral 

magnet without disorder 

oc N2tt^S'^ (l + (J(q + q' _ t) + J(Q _ q' _ r )), (44) 

Bragg T 

where N is the number of sites in the lattice, q' = t are reciprocal lattice vectors, and delta-functions describe 
the magnetic Bragg peaks (satellites) at Q = ±q' -I- r. Terms in Eq. (1551) linear in p give zero after averaging over 
disorder configurations. Most of the second order in p terms give either zero or contributions proportional to Eq. dm) 
with a small factor cQ^. The only important quadratic in p term has the following structure: 

•S' + Qz) X! VmPjm cos (to - n)q', (45) 

in,jm 



where the line denotes the averaging over disorder configurations. This averaging can be easily carried out using the 
following expression for the Fourier transform of the field from a single dipole: 


dre 


2k-] 


1 

47rJl ,k.R„ / ^/2 -ikj2\ 


I-Si 

1 

o 

1 

ft) 


f - Ro -b ^Oz 

J 

pjo " V ) 


(46) 


where Rq specifies the dipole center, k = {kx\/Ji/Jo,kyyJJi/Jo,kz), and k = {kx,ky,kz) (cf. Eq. (ITOl) '). As a result 
one obtains for the elastic cross section 


da f da \ 


— 


dVl \dQ. J 


NcS 


(> + a.) E 




1 - cos {Qz +q' -Tz) , 1 - COS {Qz - q' - Tz) 


V 


(Q + q' 


(Q - q' - T 


(47) 


where the first term is given by Eq. (1441) and the second one describes the diffuse magnetic scattering due to the 
disorder. Quite unexpectedly for diffuse scattering caused by disorder^ii the second term in Eq. dm has the power- 
law singularities at the magnetic Bragg peaks positions (cf. Eq. (1551) 1. Then, one obtains that Bragg peaks acquire 
power-law decaying tails (see Fig. a. This feature is attributed to the long-range character of the perturbation made 
by defect bonds. 


D. Magnon spectrum renormalization in the layered magnet with DMI 

In this section, we discuss defects impact on the magnon spectrum in the model ([1]) . We remind first the well known 
results for the pure system. 


1. Spectrum of the pure system 


One obtains for the bilinear part of the Hamiltonian using Eqs. o and m 




k ■- 


Akajak - ^ (okO-k + 


.D^ 


Ak = 25'(Jo(l — cos/cz) -I- Jl(2 — cosfca; — cosfcy)) -I-S'—;-(2 — cosfcz), 

2Jq 


(48) 

(49) 
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FIG. 4: Sketch illustrating Eqs. (1471) and (1851) for elastic neutron scattering cross-section at systems with defect bonds. The 
magnetic Bragg peak at momentum transfer Q = q' is shown by solid line (Q^: = Qy = 0). Power-law decaying tails are shown 
by dashed lines which are given by the second terms in Eqs. (021) and (IMll . 

(a) 

■ 


FIG. 5: Diagrams giving leading corrections to the magnon spectrum in the first order in the defects concentration. 
Then, the bare gapless spectrum reads at small k as 

4 °’ = + Ji(fc 2 + fc 2 )][Jofc 2 + Ji(fc 2 + ^ 2 ) + ( 5 I) 

Two regimes can be distinguished 

4°^ = SDk, k<^D/Jo, (52) 

4°^ = SJo'k^, D/Jo < < I, (53) 

where k = ^k^J- (fc^ -I- ky)Ji/J q- 

8. Spectrum corrections 

We imply first that only DMI is changed at imperfect bonds. One has to take into account diagrams shown in 
Fig. [ 5 ] to find the spectrum corrections. Calculations are simplified by the fact that vertexes in all of the diagrams 
are proportional to Udm which is much smaller than exchange constants. Some details of the cumbersome diagram 
analysis can be found in Appendix where the following expression for the magnon energy is obtained: 

fck = Scqudm.{2 - coskz) + c-^(/ik + hk + hk), (54) 

ao 

where /ik, l 2 k, and / 3 k are smooth functions of k which are of the order of unity and which are given by Eqs. (IA14I) . 
(IA19|) . and (IA23I) . respectively. 

The magnon damping is given by the following term which stems from the diagram shown in Fig. [SJb): 

{SudmDY t 

7k = c 




Ek JqJi 27r 


(55) 
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where t = 1 and 1/2 for, respectively, k <C D/Jq and k D/Jq (see Eqs. 

It is seen from Eqs. (ESI) and (l?T]) - (l551) that the magnon damping is much smaller than the bare spectrum at all 
momenta. In contrast, the correction to the magnon energy (15411 (that is finite at A: = 0) becomes much larger than 
the bare gapless spectrum at sufficiently small k. It implies that our results obtained in the first order in c are 
inapplicable and further analysis is required for particularly small k that is out of the scope of the present paper. 
Besides, this effect can be screened in real materials by a small gap in the spectrum originating from, e.g., a small 
anisotropic interaction. 

Let us take into account also the defect in the exchange interaction on imperfect bonds. In contrast to the defect 
in DMI, one cannot assume in general that \ufm\ Jo-, Ji- Then, one has to sum an infinite set of diagrams of the 
type shown in Fig. [SKd) to find spectrum corrections in the first order in c. As a result of tedious calculations some 
details of which can be found in Appendix lAl one leads to the following counterparts of Eqs. (IMl) and (1551) : 


fck = Scq (udm - (2 - cos fcz) + c 

qUex\'^ t 

7k - c 5 [udm 2 J JnJ, 27 r’ 


Q'^ex) j 2Q 

- j -Ak + C-—/2k + C 


2 { Udm - qU ex )Q 

Jo 


hk, 


where Q is given by Eq. 


(56) 

(57) 


E. Some other layered helimagnets with DMI 


In this subsection, we discuss briefly some other models of layered helimagnets with defect bonds to which our 
theory is applicable after some modifications. The first model differs from that discussed above by the sign of Ji 
(i.e., by replacement in Eq. o of Ji > 0 by — Ji < 0) that results in AF xy planes. It is convenient to rewrite 
transformation @ in the following more general form: 

= 5S;cos(koR.„)-5)(;sin(koR.„), 

Sfn = 5f;sin(koR„)+5f,;cos(koR.„), (58) 

qz _ qz' 

where ko = (tt, tt, q) is the vector of magnetic structure. In pure system, the spiral propagates along z axis and q is 
given by Eq. ©. The operator of perturbation (l30)) and the system of equations for pin dSD, which determine the 
spin texture around defect, remain the same. Then, Eq. (1551) for the corrections to q does not change either. 

Let us assume that both Jo and Ji are antiferromagnetic (i.e., we replace in Eq. m Jo > 0 and Ji > 0 by — Jo < 0 
and — Ji <0, respectively). In this case, ko = (7r,7r,7r — q) in Eq. (l58l) . The operator of perturbation ([30|) from a 
single imperfect bond changes its sign that leads to the dipole with opposite dipolar momentum. Then, the vector of 
magnetic structure acquires the form (tt, tt, tt — g — Sq), where Sq is given by Eq. (I39|) . 

The third model contains antiferromagnetic triangular xy planes (i.e., we replace Ji > 0 by —Ji < 0 in Eq. © and 
assume that xy planes are triangular). This model is relevant to CsCuCla (see, e.g., Refsj ^^d^ and references therein). 
In pure system, 120° spin ordering is realized in each xy plane and a spiral ordering arises along z axis. Then, the 
vector of the magnetic structure ko can be equal either to (0, ^,q) or to (0, —^,q) (we assume for simplicity that 
distances between all nearest neighbor spins are equal to unity) which describe 120° spin structures with different 
arrangements of chiralities of triangles in xy planes (see, e.g.. Ref.—). One obtains the same operator of perturbation 
dSOD. The system of linear equations has the form (EU, where J\ should be replaced by Ji/2 and one has to take into 
account that there are six nearest neighbor spins in xy plane. In the continuum limit, we obtain Eq. (1181) in which Ji 
should be replaced by 3Ji/4. Counterparts of Eqs. (1381) and (l39l) have the form 


Q 


6q 



i^dm '^ex 



3Jo “1“ 3(/i + '^ex 


Jo + Jl 


O I Jn 4_ 2. ( / 3 _ / 3Ji A 

371-Ji it VV V J 


C 1 Udm Uex 


D 

To 


To + Jl + 


1 + Jl/Jo 

Jq \ f j _ / 3 Ji ^ 

SttJi ' TT \ \J 4Jo y 3Ji+4Jo j 


< 1 . 


(59) 

(60) 


The vector of magnetic structure has the form (0, ±^, g + 5q). Changing of the sign of Jo in this model leads to the 
replacement of Q by — Q and to the vector of magnetic structure (0, ±^, tt — g — 5q). 
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Results for the spectrum corrections in these systems are not simple modifications of those obtained above for the 
ferromagnetic exchange because all these models have different bare spectra. Corresponding calculations are out of 
the scope of the present paper. 


III. CUBIC B20 MAGNETS 
A. Pure cubic B20 magnets 

Our consideration of cubic B20 magnets is based on Refs j^^’^° which are devoted to pure systems. We present in 
this subsection the well-known results which are important for further analysis of disordered systems. For discussion 
of low-energy dynamics, the following Hamiltonian is proposed which contains the exchange coupling T-Lex, DM term 
TLdm, and small anisotropic exchange interaction (AEI) Hae- 

"Ho = 'Hex +'Hdm +'Hae, 

'Hdm = — ^ DrR' • [Sr X Sr>], 

where summations on R and R' are taken over all sites of a simple cubic lattice and u = x,y^ z. As it is frequently 
done in theoretical considerations, we take the cubic lattice structure rather than the full B20 structure mainly for 
technical simplicity. Besides, a little is known now about interaction between four magnetic ions in the cubic unit 
cell of the widely discussed itinerant materials having B20 structure. It is assumed that all interactions in Eq. m 
act between nearest neighbor spins: Jrr' = J, Drr/ = D, Drr<||(R — R'), and Err/ = F. We imply below that 
J ^ D ^ F and put the lattice constant to be equal to unity. The following local orthogonal coordinate frame is 
defined at each site: 

Cr = a cos(q • R) -I- b sin(q • R), (65) 

T 7 r = b cos(q • R) — asin(q • R), (66) 

Ir = c, (67) 

where o x b = c. Spins are represented in the local coordinate system as Sr = -S^Cr + + *5'r^r. We use the 

Holstein-Primakoff representations (l5|) for spins components with the following axes correspondence: x' O 

y' O rj, and 0 O 

The ground state energy per unit cell has the form at g <C 1 

Eel = -JS^ (^ ■ y) ■ 

where I = + b^). Obviously, Ed is minimal if q||c, i.e., spins rotate in the plane perpendicular to q. The 

direction of q relative to the lattice is determined by the last term in Eq. (1551) . Eor E > 0, q should be directed along 
the cube edge to minimize the cubic invariant /. If E < 0, one infers that the helix vector is oriented along one of the 
cubic space diagonals and I = 2q^j3. In both cases one has 

D , , 

q=yc. (69) 

The main role of AEI is to determine the q direction and it can be omitted in other calculations due to its smallness. 
As E < 0 in many B20 magnets including MnSi, we discuss this case below. Henceforth, c is directed along one of 
the cubic space diagonals. 

The bosonic analog of spin Hamiltonian dSl) has no terms linear in Bose-operators and one has for the bilinear 


(61) 

(62) 

(63) 

(64) 
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terms 



n 


( 2 ) 

dm 


JS (a^flR + a^^g_^aR+eJ f 1 

R,i/ L \ 


(®R®R.+e^ + flRaR+e^) ^ 

[(a^aR + a^+g^CR+eJ - 

R,i^ 



+ (oRCtR+e^ + ®R®R+e^) 


2 (oROR+e^ + ®R®R+e„ + OrOR+g^ + ORO^^g^) , 


(70) 

(71) 


where v = x,y, z and are basis vectors of the cubic lattice. 


B. Perturbation of the magnetic ordering by defects 

Let us consider an imperfect bond between sites Rq = (0,0,0) and Ri = Gz = (0,0,1). The perturbation in the 
Hamiltonian has the following form: 


T — ^dm ~t“ Vea; — ^dm (®z * [®Ro ^ ^Ri]) ^ex^Rg ■ Sr^ . 

Omitting terms containing products of more than two Bose operators, one derive for Vdm and Ve: 


(72) 




Tea, — SUq 


ttg ao + a^ai — - (ooai + Oq a^*" + Og Oi + Oq) 


(oo^oo + a)^ai) ( 1 ~ ^ ) + (nooi + 00^01")% ~ (ajoi + a('’ao) ( 1 ~ ^ ) + 


Su. 


dm 


V3 


(oo + ad" - Oi - a(^), (73) 


(74) 


+ Qz\l -^(oi + aj" - Oo - aj) 


where indexes 0 and 1 stand for Rq and Ri, respectively. To dispose of terms in the Hamiltonian linear in Hose- 
operators, we make the shift similar to (nm which we write in the form 


Or — ^R + Pr — 6r -f Pr + ip'^, 


(75) 


where p^ and p^ are real. Simple but tedious calculations show that the following conditions should hold in order 
terms in the Hamiltonian vanish which are linear in operators 6r and b^: 

E 


j ( 2pR - PR-e,. - PR+e,. - q„p-R + if (PR+e„ + PR-g,.) + q^9(2pR “ Pk+G„ “ P'll-eJ 


{qzUex — Udm/\/3)(5R,Ro — (5r,rJ — 


ql , ~ 1 

Uex[p-R.o - PRi + y(PRi - PRo)) + -^UdmqipRo “ PRi ) 


JR.Ro 


ql , ~ 1 

Mea;(PRi “ PRq + y (PRq “ PrJ) + -^UdmqiPRi “ PRq ) 


3R,Ri 


VR. 


(76) 


Imaginary parts of equations (1751) form a linear homogeneous system of equations for p^ which gives Pr = 0. Real 
parts of Eqs. (1751) give equations for p^ which have the form similar to that of Eq. (1571) 


E \ 3^®) ~ Pr—G„ ~ Pr+g„) 


= (<5r,Ro ~ <Jr,Ri) 


{qzUex l^dm/V^) “f ^^ex 2 ^ ^'^dmq^ (PRq PRi ) 


Eqs. (1771) can be solved as Eqs. (15T1) with the result (cf. Eq. (1^ 1 

1 


p'(r) = J|E 


2 47r \ |R — Ri 


|R — Ro 


VR. 


(77) 


(78) 
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FIG. 6: Defects in cubic magnets considered in the text, (a) Imperfect bonds (marked by dashed lines) which can be oriented 
along each cubic axis with equal probability, (b) Defect with six imperfect bonds which can be relevant to mixed compounds 
of the type Mni_2;Fe3;Ge at a: <C 1 or a: « 1. 


Then, is given by the field of a dipole which momentum has the form 


where t = 2 — {^/2 — 9/4)/7r si 2.27 and 


a 

P 



7+ tp 


'^dm 

V3J(1 + gV6)’ 

'^ex ~\~ UdmQ/'i UexQ /6 

J 


(79) 

(80) 


(81) 

(82) 


Let us turn to the system with a finite concentration c <C 1 of such defects. We assume that randomly distributed 
imperfect bonds orient randomly along three cubic axes as it is illustrated by Fig.[6Da). In this case, a finite correction 
to the spiral pitch arises. In terms of the electrostatic analogy, the system “polarization” P = cd/\/Zt arises that is 
directed along the cubic space diagonal. Correspondingly, the correction to the spiral pitch q has the form 


5q 



(83) 


A substantial reduction should be pointed out of defects impact on the system properties at Udm ~ Uex<l that follows 
from Eqs. (1501) and ((5T]) . 

To verify Eqs. (1501) -(155 )) . we perform numerical calculations for a set of model parameters. We minimize the classical 
energy of clusters with open boundary conditions containing up to 100^ sites in the following way. Starting from a trial 
configuration, we arrange all magnetic moments along their current molecular fields. After performing this procedure 
many times (^ 10® -F 10^), the system stabilizes and we take the Fourier transformation of the final configuration 
(ignoring spins near the cluster boundary) which has a peak at the spiral vector q (for the given disorder realization). 
Averaging over 10 -F 20 disorder realizations, one obtains the spiral vector. Representative results of such calculations 
are shown in Fig. [T) It is seen that the agreement is excellent at c < 0.03 of numerical findings with Eqs. (1501) - (155)1 . 

Another type of defects in B20 magnets which we consider is presented in Fig.[0Db). It looks more natural for mixed 
compounds Mni_a;Fea;Ge considered recently experimentally: one expects that substitution of one magnetic atom by 
another in a unit cell of itinerant material changes couplings of this unit cell with all its neighbors. The system of 
equations (HD) describing distortion of the spiral ordering caused by one imperfect bond is linear. Then, the result 
for the considered type of defect is a linear combination of solutions for six defect bonds shown in Fig. (Db). As a 
consequence, the “polarization” and the correction to the spiral pitch are six times as large as those for one imperfect 
bond: P = 2\/Zcdc and 


5q = 2\fZcQ 


( 84 ) 
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FIG. 7: (Color online.) Correction 8q to the spiral vector as a function of defect concentration c for J = 1, D — 0.3, u^x = —0.2, 
and Udm = —0.6. Analytical result is given by Eqs. (I80I) - (I83II . Numerical result is obtained as it is discussed in the text. 


(cf. Eq. (1551)). 


C. Elastic neutron scattering 


For the model of imperfect bonds shown in Fig. IHKa), the main difference from layered magnets discussed above 
is that there are dipoles with momenta directed along three cubic axes. The concentration of dipoles directed along 
each cubic axis is c/3. Taking this into account, we get the following expression for the elastic neutron scattering 
cross-section after tedious calculations (cf. Eq. (ITTl) ): 


da f 



T iy—x,y,z 


1 - COS(Q^ + q'v- Ty) 


1 - cos(Q^ - g/ 





(85) 


where the first term is given by Eq. dBl). Then, Bragg peaks in B20 magnets acquire power-law decaying tails (see 
Fig. S]). The last term in Eq. (1551) should be multiplied by 6 in the case of defects shown in Fig. (IKb). 


D. Magnon spectrum renormalization 


It is well knowi>i2i2£ that the bare magnon spectrum obtained from Eqs. dZOl) and dm) has the form 

= SJqk, k D/ J, 

4°^=S'Jfc^ D/J<^k<^l. 


( 86 ) 


There is also a small gap in the spectrum which can be a result of magnon-magnon and magneto-elastic interactions.— 
This gap is important for interpretation of some experimental data obtained in B20 magnets—^— 

Carrying out calculations similar to those for layered spiral magnets, we obtain for corrections to the spectrum 


fek 

7k 


E Scq / qUex\ 1 \ , 

[Udm -^ j (2- COS k^) + c 

‘77 ' ■ 

c —b 

V 


(Udm - qUexT J , T . i'^dm - qUex)Q j 

-3j- 


[ '^dm 


IStt’ 


(87) 

( 88 ) 


where v = x^y, z and Jik, d 2 k, and / 3 k are values of the order of unity which are smooth functions of k having the 















15 


following form: 


hk 


hk 


16 (27r)® 


uKi (IK2 ■ 


1 


•^k '^ki ^lt.2 

- I dkidk2 - l-COs(fcl^ 3 +fc 2 . 3 ) 

,3 j (ek-Eki-ek 3 )(ki+k 2)4 

(1 + COS + /c2i/J — cos— cos k2i,^){l + cos + ^21^2) ~ cos A:ii,2 — cos k2u2) 
J 1 /* ji ji ^ COs(fci^-^ + k2u^) — COsA^ii/^ — cos ^21/1 . 2 ^1^2 k2i/2 

I ^kldk2-^^-’ 


^ (27r)6 


(ek - Eki - Eks) (ki + k2) 


(89) 


(90) 

(91) 


where 1 ^ 1 , 2,3 = x, y, z. One can see from Eqs. (IMl) and (IM)l that the damping is small compared to the bare spectrum, 
7k The correction to the magnon energy (1871) is much larger than the bare spectrum (1861) for small enough 

momenta signifying a new physics at such k. However, this effect is screened in real B20 materials at c <C 1 by the 
small gap in the bare spectrum mentioned above. 


IV. SUMMARY 

To summarize, we develop a theory describing spiral magnets with bond disorder at small concentration c of defects. 
It is assumed that both DMI and exchange coupling are changed on imperfect bonds. We obtain qualitatively the 
same physical picture in two models which are considered in detail: layered and B20 cubic helimagnets. Using the 
Holstein-Primakoff spin representation, we find the distortion of the spiral magnetic ordering around a single imperfect 
bond. It is shown that values of additional turns of spins caused by the impurity are given by Poisson’s equation 
for electric dipole. Thus, the magnetic ordering distortion from a single imperfect bond is long-range: values of the 
additional turns of spins decay with the distance r to the defect as 1/r^. Poisson’s equations for the dipole in the 
corresponding models on lattices with space dimensions d>2 give the power-decaying law 

At finite concentration of randomly distributed defect bonds, we calculate the observable quantities by averaging 
over disorder configurations. We find that the direction of the spiral vector does not change and its modulus acquires 
a correction Sq given by Eqs. (l39l) and (IMl) in the two models considered. For defects of the type shown in Fig. EDb) 
in cubic magnets, 6q is given by Eq. (15^ . It is seen from these equations that the spiral vector correction can be 
zero, positive or negative depending on the particular parameters of defects. For negative 6q, the sign of chirality can 
change even at c <C 1 if defects are strong enough. 

In the elastic neutron scattering cross-section, defects manifest themselves in two ways. First, magnetic Bragg peaks 
(satellites) are shifted from reciprocal lattice vectors by ±(q-l- dq) (i.e., by values defined by the new spiral vector). 
Second, diffuse scattering arises which has power-law singularities at Bragg peaks positions. Then, each Bragg peak 
acquires the power-law decaying tails (see Eqs. gZl), dill), and Fig. n. This feature is attributed to the long-range 
character of the perturbation made by defect bonds. 

Corrections to the magnon energy and to the damping caused by scattering on defects are given by Eqs. (Unj-dlZl) 
and ((87l) - ((88)) in layered and B20 magnets, respectively. The magnon damping is found to be much smaller than the 
bare spectra in both models. Although magnons are well defined at fc ^ g in both models, the ratio yk/^k ^ c/k is 
quite unusually large. Remember, this ratio is normally proportional to a positive power of k and it does not exceed 
c in magnetically ordered gapless magnets (see, e.g., Refsi^^— and references therein). However we have obtained 
recently that 7k/£k ^ c/fc^ under certain conditions in gapped phases of 3D spin systems with bond disorder^ 

Corrections to the magnon energy exceeds the bare spectra at small enough momenta. This signifies that the 
analysis cannot be restricted by the first order in defects concentration at such k. It can also signify a localization 
of long-wavelength magnons (see, e.g., ReUSI and references therein). Consideration of this point is out of the scope 
of the present paper. Besides, these small-energy peculiarities can be screened by a small gap in the bare spectra 
originating from a small low-symmetry spin interaction. 

Although all calculations for layered helimagnets are performed for the model with EM exchange interactions, 
the results obtained (except for the spectrum renormalization) are applicable after simple modifications discussed in 
Sec. Ill El to many other layered helimagnets with bond disorder. 

Our consideration can be relevant to Mni-ajEe^Ge considered recently experimentally in Refi^. But we are unable 
now to verify our theory due to very small amount of experimental data at a: ~ 1. For instance, there are only 
three experimental points on the plot for dependence of the spiral vector modulus on a: at a: > 0.75. Then, further 
experimental activity is needed in this field. 
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Appendix A: Calculation of the magnon spectrum renormalization in layered spiral magnets 


First, we take into account the imperfection of the DMI only. In addition to terms in Vdm presented in Eq. (nni, 
one needs also the following terms for the magnon spectrum calculation: 

^dJi ~ + a'^m+l'^hm+l — 2 {'^hm^hm+l ^hmahm+1 + ahm'^hm+l + adrn+l^hm)j , (Al) 


hm 

^dm ~ [ CL^m^hmO'hm+1 + ^hm'^hm+l^hm + 

hm \ 


+ 2 +2 
^hm+l^hm+_i ^m+l^f^rn+1 


4 4 

+ ^2 „+2 




+ „ + 


O-hm+l^hmO-hm+l ^hm^hm+l'^hm+l 


_ ^hm^hm _ \ 


(A2) 


4 J ’ 

where sums run over sites involved in defect bonds. Besides, one has to take into account terms in the Hamiltonian 
containing products of four Bose operators 

E —^ { ~\~2 ~\~2 ——|— 2 —f— 2 \ 

^in^in+l^inClin+1 — -^[ain+l<^inCiin+l + ainClin+l + 0,in^in+l^in+l + ^in^in+l^in) 


{ij)n 


(A3) 


where we omit terms of the second order in D/ Jq <C 1. 

Eq. (IA3I) gives the following terms after shift (TT^ which contain products of one operator of creation and one 
operator of annihilation: 


= -2EE^. 


in j 


^in^iniPjn PinPjn) Pjn) 


(A4) 


where j enumerates nearest neighbors of the f-th site in the n-th plane, Jj = Jq and Jj = Ji for neighbors from 
different planes and from the same plane, respectively. It can be shown that terms containing products of two 
operators of creation or two operators of annihilation give a negligible correction to the spectrum. Introducing the 
Fourier transform 


bin — 


//V 


E^ke 


— ik R-i, 


(A5) 


we 


have for Eq. (IA4l) 




,2Rin-(ki-k2) 


ki,k2 


in J 


Pjn PinPjn 2 (/^m Pjn^ 


(A 6 ) 


where k 2 j = k 2 z and k 2 j = k 2 x or k 2 y for neighbors from different planes and from the same plane, respectively. 

In much the same way, one obtains for terms containing products of three Bose-operators and stemming from "^4 
(Eq. dMll) 


^ 

4 ^ 3/2 


E E E ( 2 p,„ - p„) - (l + + h.c. 


ki,k2,k3 


in J 


(A7) 

Taking into account only terms in Eq. ini containing products of one operator of creation and one operator of 
annihilation which give the main contribution to the spectrum renormalization, one obtains 


E 


ki,k2 


hm 


I _|_ ^K^l^-k2z) _ 


^ikxz _j_ g ik 2 z 


(AS) 
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We have from Eq. 


after the Fourier transformation 




iV3/2 


ki,k2,k3 


hm 


(A9) 

Let us start the spectrum calculation with Eq. (IA8E As q\udm\ ^ >^o,ij the main corrections to the magnon energy 
i5ek and to the damping qt originate from diagrams shown in Fig. [SJa) and[Sl[b), respectively, which give 


= Scqudmi^ - coskz), 


(AlO) 


( 1 ) 

7k = S 


{SqUdm)^ 1 

A^2 ^ _ iQ 

{SqUdmf 

(27r)3 




,ikz 


d^ki 




1 ^^{kz-kiz) _ 


{SUdmDY t 

Ek JqJi 27r’ 


^—ik\ 


hm,h'm' 


•(ki-k) 


(All) 


where S denotes imaginary part, t = 1 and 1/2 for, respectively, k <C D/Jq and k :#> D/Jq (see Eqs. (ICT) " (I551) 1. 
hereafter the line over an expression denotes averaging over disorder configurations, and we take into account that 
only terms with h = h' and m = m' survive after the averaging over disorder configurations in the double sum over 
hm and h'm'. 

The main contribution to the spectrum renormalization from Eq. (IA9I) originates from the diagram presented in 
Fig.ISc). After integration over internal frequency and averaging over disorder configurations, we have for it 


cu, 


2 ^ \ _ 

dm g^2 2-j _ £ _ jQ 

ki,k2 


+e*'=^^)(l-e-*'=0 + 


gi{ki^+k2^-k^) _ ^ 


(A12) 


The imaginary part of this equation is of the order of cu'^^k^^f ei^. Thus, it is larger than Eq. (IA11|) only for quite 
large momenta, k ^ {DfV To The correction to the magnon energy has the form 


Se 


( 2 ) 


7ik = 


„ ^dm j 

C^—71k, 

7o 

s Jo 


dkidk: 


1 - cos(fci 2 ; + k 2 z) 

£^k £^ki ^k2 


16 (27r)6 

that should be taken into account together with Eq. (lAlOp . 

The main correction to the magnon energy from Eq. (IA6I) is given by the diagram shown in Fig. [SKa): 


(A13) 

(A14) 


(cos fcj l)(pm PinY'- 

in j 


(A15) 


It is negligible compared to Eq. (jA13l) being of the order of cu^^k"^. Contribution to the damping from Eq. (jA6l) 
stems from the diagram depicted in Fig. [SKb) and it has the form 


"MAE 


^ £k- £ki -iO . 

ki in I'n 


Jj Jj/eFR-n R,'n')'(ki k) I^p2^ _ p2^-^ (^p2^^ _ p2^^^ (a16) 


It can be discarded being of the order of cu^^k^/e\^. 

The loop diagram shown in Fig. [SKc) with three-particle vertex (IA7|) gives 

E JJj 2 \ A (It cos (^k\j^ “t“ k^ji ) — cos kijj^ — cos k 2 ji ) (1 T cos (/^ij 2 T ^ 2 j 2 ) — cos k^j^ — cos k 2 j 2 ) 
AAf3 2^ 


4iV3 

71J2 ki,k2 


£k £k.i £k.2 7 O 


E 


^ y PinPjm^ 
in,jm 


i(Rin—Rjm)-(l^l+k2) 


(A17) 
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Using Eq. (l46l) . one finds the following valuable contribution to the magnon energy from Eq. (IA17I) : 

,2 




(A18) 


~ 2^ 4Jo (27r)6 


dki(ik2 


J1J2 


1 - cos (fciz + k2z) 

(ek - £ki - ek2)(ki + ^2Y 


X (1 + cos + fc 2 ji) — cosfciji — COSfc 2 ji)(l + cos (fcij 2 + feja ) “ COS — COSfejs), (A19) 

which is of the order of as Eq. (|A13E The imaginary part of Eq. (IA17I) is of the order of /sk. Thus, it 

is larger than Eq. (lAllI) only for quite large momenta, k {D/. 

There are also corrections from the diagram shown in Eig. [5l^b) which is built using both and . The 
corresponding expression has the form 

E JyqUdmS _ I /jRi„.(ki-k) 

/V2 2^ ejj. — — jO 


x^e 

jm 


iRjm'Ck—ki) 


{in} 


I _|_ gi(kiz-k^) _ 




'2 _|_ 
Pjm+e„ ~ PjmPjm+e„ ~ Pjm—e^Pjm ~r C 


—ik\u 


PjmPjm-\-ei, 


Pjm+eu Pjr) 


+e 


ikiu 


PjmPjm—eu 


P^ + 

rjm—ei^ ' rjm 


Th.c.y 


(A20) 


Corrections to the magnon energy and to the damping from this expression are negligible being of the order of cu^^D 
and cu^^Dk^/sk, respectively. 

The second correction of this type comes from the loop diagram presented in Fig.jSJc) which contains both vertexes 


and Hf and has the form 


,(3) 


dm 


E 


2 ^dmJu 


^ a; - Ekr - Ek-, + ^0 
«1,«2 




[S' 


,iRin-(k—ki —k2) 


gi(kz-kiz-k2z) _ 2 


x^e 

jm 

H” Pjm-\-ei^ 


2Rin-(ki+k2-k) j 


Pjm 


^—iki, 


:,ikiu I f,ik2u 


-iku 


-(l + e 


i{kx„+k2v)'\ 


(A21) 


^ — iku 


1 - 


e^kiy + gife2. N I 




+ h.c. 


The imaginary part of this expression at small k is of the order of cu^^k^/ei^. Then, it is negligibly small. The real 
part is given by the following equation: 




Ji 


hk = 


SJ, 

2(27r)6 


1 + cos(A:i^ H- k 2 u) ~ cos kh, — cos k 2 iy . 2 
akiak 2 -^- — -sin 


(ek - eki - ek2) (ki + ^ 2 ) 


(A22) 

(A23) 


which is of the order of and should be taken into account. 

Let us take into account the defect in the exchange interaction (ED which has the following form after the Fourier 
transformation: 


V(2) = Sue, E E (1 - 6*'=-) (1 - 6-*'=-) , 


(A24) 


{m} ki.k2 


where {in} denotes imperfect bonds. In general, one cannot assume that |uea;| “C Jo,Ji as it was for Udm- Then, one 
has to sum an infinite set of diagrams of the type shown in Fig. [SKd) to find spectrum corrections in the first order in 
c from Eq. (IA24I) . As a result, the Green’s function denominator has the form 


G(a;,k)-i = cc-4°)-T(w,k), 


r^k) - c2S„„(l cost.)(l 2S~=/ 


dq (1-cosq^) 
(2’’")^ oj — — iO^ 


(A25) 

(A26) 
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One has from this expressions 




(A27) 

7k 

+ huexY ’ 




fck 

1 “1“ ^4^ex 

(A28) 

h 

oo /" 1-cosg^ 

J (27r)3 AO) • 

(A29) 


These results are negligible compared with those stemming from the defect in DMI which are considered above. There 
are also corrections from diagrams of the type Fig. [SJd) made both from Eqs. (IA8I) and (IA24I) . Their analysis shows 
that they are also small. 
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